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Abstract 

It has been shown that it is possible to extract values for critical couplings and ^string in matrix 
models by deriving a renormalization group equation for the variation of the of the free energy as 
the size N of the matrices in the theory is varied. In this paper we derive a "renormalization group 
equation" for the Penner model by direct differentiation of the partition function and show that it 
reproduces the correct values of the critical coupling and 7 string and is consistent with the logarithmic 
corrections present for g = 0,1. 
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It was observed by Brezin and Zinn-Justin Q that the known double-scaling behaviour of the free 
energy F of the matrix models associated with c < 1 theories 

Fa A 2 - 70 /(AiV 2/71 ) (1) 



where A = g — g c , 70 + 71 = 2 and 71 = (1/12) (25 — c + — c)(25 — c)) meant that they should satisfy 
a renormalization group equation of the form 

N4r T -m4-+7(9))F = r(9) (2) 
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where 7(5) = 2, with the critical coupling(s) determined by f3(g c ) = and the scaling exponents by 
71 = 2/(3'(g c ). Such an approach, considered before in other contexts (2), is of interest as a possible 
means of extending matrix model calculations beyond the so-called c = 1 barrier. They showed that a 
lowest order calculation, in which a row and a column were explicitly integrated out in a matrix model 
at finite N, gave values that were reasonable and suggested that pushing the calculation further would 
improve the numerical accuracy. However, a direct continuation of this approach proved disappointing, 
with no sign of a convergence to the known critical couplings Similarly, applying the methods to a 
d = 1 matrix model gave results that were correct qualitatively but not very accurate quantitatively [|J . 
The vital missing insight came from work on the renormalization group for vector models ||, where it 
was observed that it was necessary to take account of the reparameterization invariance of the theory via 
the Schwinger-Dyson (ie loop) equations in order to eliminate unwanted induced couplings. The same 
approach gave excellent results when applied to d = matrix models as well || , whether in a eigenvalue 
representation or by direct integration of a row and column in the style of Brezin and Zinn-Justin. 

If we integrate out a single eigenvalue A in the style of j(| to take us from a JV + 1 x JV + 1 matrix 
model to a N x N matrix model then we hnd, symbolically, a modified action 

Z N ~ J d\ J d Ar2 $exp[-(7V + l)iry($)-(iV + l)y(A) + 2<rlog|A-$|] (3) 

where V($>) is the original action. This in turn leads to a renormalization group equation as N — > 00 for 
the free energy F 

N M + 2 ) f = ~n~ < trV{q>) > +V{ ~ X) ~Jf <tr log |X " $l > + (4) 

where A is the saddle point value of A arising from carrying out the integration in equ(^). The use of the 
loop equations, which encode the reparameterization invariance of the theory, enables one to eliminate 
the complicated induced interactions arising from the logarithmic term and get back to a beta function 
that depends only on the original couplings 
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at the expense of including terms nonlinear in dF /dg in G. The difficulties arising in the application of a 
renormalization group procedure to the d = matrix model reside essentially in the induced logarithmic 
term that is produced by the integration over an eigenvalue (or, equally, a row and a column). It is thus 
not unreasonable to ask what happens in the case where such a term is already present in the original 
action. Such models have been considered already as matrix models for open strings, as the logarithmic 
term has the effect of generating boundaries in the worldsheet fjj. In these open string models the 
couplings in front of the logarithm and in its argument are fed in by hand, rather than being determined 
by the original closed string action. 

The simplest non-trivial action that can be concocted with a logarithmic term, namely a linear plus 
logarithmic action, has been studied quite extensively in the context of calculating the virtual Euler 
characteristics for the moduli space of Riemann surfaces S as well as for being an interesting specimen 
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of a matrix model from the string theory point of view @, [T(| . If we use the notation of Distler and Vafa 
[Tcfl , the partition function we wish to consider is 

Z = f d N2 $expNttr[log(l -$) + $] (6) 

where the integration is taken over positive definite M = 1 — $ to avoid problems arising from the branch 
cut in the logarithm. Penncr showed that the free energy F — log Z for this model could be expanded as 

f = Y J N 2 - 2 n 2 -^~ n x 9 , l (7) 
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where the coefficients Xg,m the virtual Euler characteristics of the moduli space of Riemann surfaces of 
genus g with n punctures, could be calculated explicitly for all g, n 

_ (-l)*(2 g -3 + n)l(2 g -l) 
Xff '" " {2gjH\ B29 (8) 

with the Big being the Bernoulli numbers. Thus we know the expansion of the free energy for all g, n in 
this case. 

We can make use of this fact to derive a "renormalization group equation" for the Penner model by 
simply differentiating equ.(^) directly, instead of having to integrate out eigenvalues or rows and columns 
in a matrix. If we reintroduce the customary factor of 1/iV 2 in front of F and write 

oo 

F 9 = E N ~ 29t2 ~ 29 - n X 9 ,n (9) 

71=1 

then it is a simple matter to show that 

N lw- t ^m + 2F ° = £ nN - 29 * 2 - 29 - n x 9 ,n (io) 

n—l 

where we have written explicitly the limits on the n summation. If we do not divide by N 2 the 2F g drops 
out. The right hand side of equ.(fl0|) is almost equal to dF g /dt which can be written, after a shift in the 
summation variable n — > n + 1 , 

dF g ^ v - 2gf 2-2 9 -n (-1)"(2.9 - 3 + n)\{2g - 1) 
-m=h 1 (2 fl )l(n-l)l 

Only the n = 1 term in the summation is missing, so if we put this back by hand we have, with the 
proviso g > 1, 

£ N- 2 H 2 - 2 v- n n Xg , n = - N- 2 ^- 2 ^ (12) 

n=i y 

SO 

N M {t + ^ + 2F * = -N- 2g t'- 2 ^ (13) 

which is a renormalization group equation of the form in 0. The critical coupling t c of the model is 
determined by the zero of the beta function, j3{t c ) — 0. As 0{t) = t + 1 by inspection of equ. ([l3|) , this 
gives t c = —1, which is the correct result. Similarly, we should have 

^ ^ WS (14) 

which trivially gives 71 = 2, thus recovering the known fact that the Penner model is a c = 1 matrix 
model. It is not surprising that we have obtained the exact values for t c and 71, as we started with an 
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exact expression for the free energy, what is remarkable is that we have not had to include any terms 
that are nonlinear in dF g /dt to write down the renormalization group equation, in direct contrast to the 
matrix models without logarithmic terms that were considered in |l . In this respect the behaviour of the 
Penner model is more like that of the vector models considered in J5|, where the renormalization group 
equation is linear, even after using reparameterization invariance to eliminate unwanted couplings. 

So far, so good but we have not yet considered the g = 0, 1 cases and the presence of logarithmic 
corrections for these. Generically, we would expect to find a double zero in the beta function when 
logarithmic corrections are present. We are safe for g > 1 as we know such corrections are absent, so 
the renormalization group equation we have written down is at least consistent. It turns out that the 
corrections are lurking in the lower limits of the n summations we have carried out to evaluate F g . In order 
to make sense of the sums for g = 0, 1 Distler and Vafa regularized by considering the three punctured 
sphere to kill the SL(2, C) invariance and the once punctured torus to kill the translational invariance 
and then integrating to get the free energies without punctures. In the matrix model we differentiate 
with respect to the renormalized cosmological constant fx, defined by 

+ £ (15) 

to generate punctures. If we reinstate the factor of iV 2 in F to compare directly with the results of Jl0| 
and differentiate the expression for Fq three times w.r.t. fj,, we find 

oo 

F (3) = -1/7V]T(1 -^/N)-^ +1 \ (16) 

n=l 

where the superscript denotes the number of derivatives. The factor of 1/N comes from the TV 2 in the 
definition of Fo multiplying 1/N 3 from the differentiations. Performing the sum gives 1 //it as N — ► oo and 
the three integrations to get Fq from Fq generate the logarithmic corrections. Similarly, differentiating 
F\ once gives 

71=1 

which sums to — 1/12/x as N — » oo. A single integration to get back to F\ again produces the requisite 
logarithm. 

The renormalization group equation is phrased in terms of the bare coupling t rather than fj,, but we 
can still carry out the summations in i 7 ^ 3 -* and F^ in terms of t. Taking the torus first we find 

^' = ^aVo < 18 > 

where the superscripts now denote differentiations w.r.t. t. By direct differentiation this satisfies 

If we differentiate our renormalization group equation w.r.t. t (remembering there is no 2F g term, and 
that the factor of N on the right hand side is now N 2 ~ 2g as we have retained the overall factor of N 2 in 
F) we get the same equation, so the results are consistent. 
For the three punctured sphere we have 

F (3) = -N 2 j- 1 . (20) 
t(l+t) 

which, again by direct differentiation, satisfies 
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If we differentiate our renormalization group equation three times w.r.t. t, use NdFjf'/dN = 2F^> and 
set g = on the right hand side after performing the differentiations (as this is ill-defined otherwise) we get 

(3) 

precisely this equation. Integrating Fq once also gives the correct value for the genus zero susceptibility 
log(l-H) — log(t). We have thus shown that, although there is no double zero for (3{t), the regularizations 
necessary for the sphere and the torus generate the logarithmic terms in the free energy and the regulated 
free energies F^\f[ 1] satisfy differentiated versions of the renormalization group equation. 

We have in a sense cheated in this paper - because the form of the free energy is exactly known we 
have not had to do any renormalization in order to get our "renormalization group equation" . If we look 
schematically at a proper renormalization group calculation in the style of where a matrix eigenvalue is 
integrated out, it is possible to see that all the new induced terms are of the form of the original action up 
to shifts in the coupling, so there is no need to introduce new couplings and then use rcparamctcrization 
invariance to gobble them up again. It is rather remarkable that a complicated potential with interactions 
of arbitrary order should give such a simple beta function and not induce any further couplings under 
renormalization. This may be because all of the possible vertices are already included in the theory. 
There could be some hints for attempts at finding a geometrical interpretation for the matrix model 
renormalization group equation for more general actions in all of this. We have found for the Penner 
model that a change in N induces a flow in the coupling constant for punctures in the surface, which 
is the only coupling in the theory for this model. The generation of an induced logarithmic term for 
other matrix models suggests that the renormalization group flow of puncture couplings might play an 
important role here too, as put forward in the speculations at the end of jy], which considered the baby 
universe structure of 2d quantum gravity. The importance of punctures was also stressed in the heuristic 
arguments of |Q , which suggested that the c = 1 barrier was the result of a condensation of "spikes" (ie 
punctures) in the Liouville theory. It should be noted, however, that the change of scale in the matrix 
model renormalization group is in the space of the matrix indices Jl3| , so a direct correspondence with 
real space renormalization and finite size scaling ideas is not obvious. 

D. Johnston was supported at Orsay by an EEC Human Capital and Mobility fellowship and an 
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model and the matrix model renormalization group. 
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